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In a continuation of ongoing Naval Postgraduate School efforts to study the 
performance characteristics of waverider-configured vehicles, a computational fluid 



dynamic (CFD) analysis of the Mach 6 




optimized Price waverider was 

Jmax 



conducted. This analysis was performed to determine the theoretical force and moment 
data over a broad Mach number spectrum, and to compare theoretical and experimental 
results in the subsonic flight regime. The CFD determination of force and moment data 
represents a continuation of the ongoing analysis of the Price waverider configuration. 
Selected viscous and inviscid flow solutions for flight conditions in the range 0.3 < NL < 
6.0 , as well as a subsonic (M, = 0.3) angle of attack sweep, were conducted using 
NASA CFD software (OVERFLOW 1 ,8b). Examination of the computed converged 
flowfield solutions suggests that the surface pressure distributions and Mach number 
contours surrounding the body are valid. Low speed force and moment coefficient data 
are shown to exhibit reasonable agreement with the available subsonic wind tunnel data. 
Additionally, supersonic CFD results show the development of the expected shock layer, 
exhibiting an attached shock bed at the design Mach number (Moo = 6.0). Evaluation of 
the computed Mach number effects on lift and drag coefficients at subsonic, transonic 
and supersonic Mach numbers suggests that the Price waverider may exhibit some flight 
instabilities across the flight Mach number spectrum. 
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I. 



INTRODUCTION 



A. BACKGROUND 



In 1 959. Nonweiler hypothesized that the pressure generated behind a planar 
shock wave might be sufficient to provide the vehicle lift. His proposal of a caret-shaped 
vehicle was the first of a class of vehicles designed for supersonic flight. The principle of 
lift generated by the pressure rise through a shock wave lead to the term “waverider”. 

The decades that followed Nonweiler’s hypothesis have brought a series of 
improvements and refinements. The latest group of proposed vehicles is a selection of 
conical flow waveriders, designed to capture and ride upon conical shaped shocks. A 
number of studies on the optimization of waverider-configured vehicles are ongoing 
[Refs. 1 and 2]. 

The Naval Postgraduate School (NPS) is pursuing in an ongoing conceptual 
development of hypersonic, carrier-based waverider-configured aircraft. Since 1992 
both individual thesis students and design team efforts have contributed to the body of 
knowledge in this emerging field of aerodynamics [Ref. 3]. Included within this study is 
the development and investigation of a proposed carrier-based interceptor aircraft, and 
among the proposed designs is the Price waverider. 

Using the NASA Ames Research Center (ARC) Waverider Code, David Price 
developed a conceptual configuration for a hypersonic, carrier-based, deck-launched. 



waverider strike aircraft, optimized for Mach=6, 




, flight [Ref. 4], The 



waverider code is a subset of the Hypersonic Aircraft Vehicle Optimization Code 
(HAVOC), created by the Systems Analysis Branch at NASA ARC. Optimization of the 
product of L/D and the specific impulse produced an integrated engine and air vehicle. 
This method is in concert with one of the fundamental design aspects of hypersonic 
vehicles. The first and most important aspect is to create and capture a conical shock. A 
second important aspect, however, capitalizes on shock compression, forward of the 
engine intake. This shock pre-compresses the engine inlet air. Aft body shaping is used 
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to expand the exhaust in an external nozzle. Thus engine air compression and expansion 
takes place external to the engine modules themselves allowing for a reduction in weight, 
complexity and cooling requirements. In this sense the engine and aircraft configuration 
must be an integrated effort, both being interdependent in all aspects of flow analysis - 
hence the need to optimize this type of vehicle using both aerodynamic parameters (L/D) 
and engine parameters (Isp). 

Follow-on work to Price’s initial development are ongoing. Based on this 
configuration, 15-inch and 8-inch root chord aluminum models were produced by NASA 
Ames for NPS research to facilitate the further characterization of the performance of the 
Price waverider. Using the smaller aluminum model, Lowell M. Johnson performed a 
series of flow visualization tests of the Price waverider in the NPS water tunnel [Ref. 5]. 
Subsonic wind tunnel testing of the Price waverider has been performed by Mark E. 
Cedrun [Ref. 6] and Michael Huff [Ref. 7], using the larger of the aluminum models. 
Tests to date indicate that flight is feasible for waverider-configured aircraft. The Price 
waverider database, while including water and wind tunnel data, lacks a Computational 
Fluid Dynamic (CFD) component. Additionally, no supersonic or hypersonic data were 
available for the Price waverider prior to this CFD analysis. 

The primary thrust of this thesis was the development of the appropriate surface 
grid for the Price waverider and the calculation of aerodynamic characteristics in 
subsonic and supersonic flight. The subsonic CFD flow solutions were aimed at 
complimenting existing wind tunnel and water tunnel data. CFD results are compared 
with the force and moment data obtained in recent wind tunnel tests for a range of 
Reynolds numbers and angles of attack. The range of subsonic, transonic and supersonic 
solutions represents a continuation of the ongoing research effort supporting the complete 
performance analysis of the Price waverider. 

B. OBJECTIVE 

The objectives of this thesis are to determine the theoretical force and moment 
data over a broad Mach number spectrum, and to compare theoretical and experimental 
results in the subsonic flight regime. The results, across the whole flight spectrum of 
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Mach numbers, are to serve as a basis for continued examination and feasibility analysis 
of the Price waver ider. 

C. PRICE WAVERIDER 

The Price waverider was optimized for Mach 6.0 flight using inviscid, conical 
shock, optimization methods contained in the Ames Research Center Waverider Code. 
The parabolic leading edge droops from the nose to the trailing edge. While the upper 
surface remains entirely in line with the freestream, the lower surface slopes to create a 
wedged-body over which the shocks are formed. Along the centerline, a pair of ramps 
precedes the engine inlet. Off the centerline the lower surface creates a 5.2 degree semi- 
vertex angle. Along the aft-body the lower surface is faired upward to meet the trailing 
edge. Waveriders derived from strictly conical flow will have a flat base at the trailing 
edge, creating enormous base drag, thus the lower surface must be faired to narrow the 
trailing edge to decrease the base drag or pressure drag of the vehicle [Ref. 8]. The fairing 
of the aft body under-surface up to the trailing edge creates a negative camber in the 
body. This negative camber is more pronounced on the centerline, where the engine 
cowling ramps back up to the trailing edge. 



3 



THIS PAGE INTENTIONALLY LEFT BLANK 



4 



n. DESCRIPTION OF CFD 



Computational Fluid Dynamic methods were developed as a means of solving the 
complex system of partial differential equations that describe a given flowfield. The use 
of computers to solve discretized equations enables the solution to be found for those 
partial differential equations for which no exact solution exists. Since 1975, numerous 
advancements have been made to refine the application of several finite element methods 
for computing flowfields; improving algorithms to improve computational efficiency and 
further defining and extending the regimes in which CFD methods can be applied with 
confidence [Ref. 3]. 

The application of CFD techniques to a known geometry is a multi-step process. 
Several software applications are used for grid generation and refinement, and flow 
solution and analysis. NPS has a suite of computer software that enables grid generation 
and manipulation, flow solving and graphic presentation of the solution. 

A. GRIDGEN OVERVIEW 

GRIDGEN Version 9 is a Computer Science Corporation developed software 
tool, which was originally sponsored for development by the U. S. Air Force. It is an 
interactive code used to create 2-, or 3-dimensional grids. The domains can be 
represented graphically within GRIDGEN or exported as data files for additional 
operations, such as providing the computational space for flow solutions. GRIDGEN 
requires the import of surface points from a previous design effort or the creation of 
surfaces manually within the bounds of the program [Ref. 9]. 

B. GRIDGEN3D OVERVIEW 

GRIDGEN3D is a batch code used in conjunction with GRIDGEN Version 9. 
This code is used to refine the grid volume by adjusting grid points internal to the block 
of points. In using of GRIDGEN3D, caution should be applied to ensure that the body 
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points are fixed to prevent the blending of body features by designating body surface as 
fixed boundaries. 

C. GRIDED 

GRIDED is an interactive, menu-driven file FORTRAN-based software 
application used to manipulate the grid. GRIDED allows the addition and removal of 
grid points or planes, the selection of indices and the performance of coordinate 
transforms. 

D. O VERFLO W1.8B OVERVIEW 

OVERFLOW 1.8b is a flow solution code developed by NASA Ames Research 
Center and is a combination of the older F3D and Chimera codes. It uses an implicit, 
finite-differencing scheme to solve the Reynolds-Averaged Navier-Stokes equations in 
strong conservative form [Ref. 10]. The input file from the user enables the selection of: 

• Flight regime or flow conditions, such as Mach number, angle of attack, yaw 
angle, ratio of specific heats, viscosity and freestream temperature. 

• Turbulence model selection (Options include Baldwin-Lomax, Baldwin- 
Barth, k-G5 and k-e models) 

• Computational controls such as time step, differencing scheme, artificial 
viscosity selection and smoothing parameters. 

• Boundary conditions (options include viscous or inviscid walls, inflow or 
outflow, or flow-through cuts). 

OVERFLOW 1.8b solutions are run on the Cray Y-MP J 94 computer, owned and 
operated by the Naval Postgraduate School. 

E. FAST 

Flow Analysis Software Toolkit (FAST) was used to compute and graphically 
represent the grid and its corresponding flow solution. The images of the waverider body 
and the Mach number and pressure profiles contained within Chapter IV are products of 
FAST. [Ref. 11] 
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III. GRID GENERATION AND FILE INPUT 



Development of an appropriate grid and the correspondingly appropriate 
boundary conditions and flowfield constraints was essential to the production of realistic 
solutions for the waverider flowfield. Additionally, a number of options exist within 
OVERFLOW1.8b, the flow solver, for refining the techniques used to optimize the 
computational efficiency or to enhance certain flow characteristics. 

A. GRID GENERATION 

Grid generation was accomplished using the previously described code, 
GRIDGEN. Using grid points, provided by NASA Ames (originally generated by David 
Price), the surface geometry of the Price Waverider was defined by importing the series 
of upper-surface and lower-surface fuselage station grid points into GRIDGEN. The 
surface domains were then built up by manually inserting the coordinates in GRIDGEN 
to create sequential fuselage station connectors. Once both the upper and lower surface 
of the waverider was completely defined by the series of linked fuselage station domains, 
these domains were consolidated into a single-domain covering the upper and lower 
surfaces. Following the initial creation of the waverider’s surface grid, a number of 
anomalies were observed in GRIDGEN in the basic geometry of the grid, so minor 
adjustments to the database provided by NASA were made. 

Those changes included a coordinate transformation, elimination of co-incident or 
nearly co-incident fuselage stations and minor body modifications to match model 
geometry. The fuselage station coordinates that were provided by NASA Ames created a 
body whose axis lay in the z-direction. Since OVERFLOW1 .8b specified M® in the x- 
direction, a coordinate transformation was performed to a right-handed coordinate system 
with the axis of the body lying on the x-axis. After an examination to ensure there was 
no loss of geometric features on the waverider surface, fuselage stations at x = 0.03014, 
0.03024, 0.03034, 0.33192, 0.63181 and 0.65947, were eliminated due to their positions, 
nearly coincident to adjacent stations. Additionally, the grid point defining the outermost 
point of the fuselage station at x = 0.93 1 64 appeared to be missing, creating a sawtooth in 
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the leading edge. A new end point was created for this fuselage station to bring the 
points at that station out to create a smooth leading edge, rather than the sawtooth 
observed in GRIDGEN. The new point was created by interpolating between the outer 
points of adjacent fuselage stations. Lastly, the engine exhaust in the wind tunnel model 
was a ramped exhaust, whereas the points provided by NASA Ames modeled the exhaust 
as a flat plate. To match the geometry of the wind tunnel model and to match a more 
realistic full-scale configuration, the fuselage station at x = 0.93816 was eliminated in its 
entirety. The resulting engine exhaust ramp geometry then mimicked the model. The 
surface grid points, with the required adjustments, are provided in Appendix A. 

The spanwise grid dimension (grid dimensions refer to the number of points in the 
grid in a given direction) of the waverider surface was set at 20 grid points each for the 
upper and lower surfaces. After some experimentation, 20 grid points in the span was 
seen as the minimum number of points required to model the surface geometry clearly. 
Additional points would increase the size of the grid and the grid density, slowing 
computations. Clustering points at the outer edge of the body’s upper surface was done 
to provide sufficient grid point density near the sharp leading edge. The body’s lower 
surface required grid point clustering at both the leading edge and at the engine cowling 
in order to ensure adequate modeling of the body geometry. The chordwise dimension of 
the grid was 32 points from nose to tail, after the removal of the fuselage stations noted 
above. Additional gridlines were added to provide higher grid density at the trailing edge 
of the engine cowling and at the trailing edge of the vehicle. 

After numerous attempts to eliminate negative cell volumes at the trailing edge, a 
cusp was added at the trailing edge to ease grid geometry. This addition extended the 
blunt trailing edge a distance of 1 .5 percent of the root chord length. This cusp is the 
only modification made that changed the computational grid shape from the model shape. 
Creation of a sharply-cusped trailing edge was believed to be more representative of final 
vehicle shape, as this addition helped to alleviate the high base-pressure drag forces 
caused by the bluntness of the Price waverider trailing edge. The final dimension of the 
grid was 45 points from nose to tail. Overall length of the vehicle, as prescribed by the 



8 



grid points provided by NASA Ames, was equal to 1.0, prior the addition of the cusped 
trailing edge. Dimensions, when quoted as a fraction of the body length, are a fraction of 
the original body length or root chord, vice a fraction of the computational model length 
or root chord. 

The resulting surface grid is shown in Figures 1 and 2. The inset of Figure 1 
shows the addition of the trailing edge cusp, increasing the total model length by 1.5 
percent. 




Figure 1. Waverider Grid Three-View and Oblique 



0.015 
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Figure 2. Waverider Surface Grid Cusp 

In addition to the complex geometry of the lower surface and the sharp curvature 
of the trailing edge cusp, the droop of the leading edge of the wing contributed to the 
difficulty in creating a single three-dimensional block that was suitable for use by 
OVERFLOW!. 8b. The vehicle configuration required a multi-step approach to creating 
a suitable 3-dimensional computational block surrounding the waverider half-body. First, 
it was necessary to split the vehicle’s upper and lower surfaces and map out the 
surrounding airspace in two separate blocks to prevent gridlines from crossing from the 
lower surface to the upper surface. Then the two blocks had to be merged to create the 
entire computational space. This was done by creating the surrounding surfaces 
including domains with the following number of grid points: 



Sting: 


20x50 


Wake Cut: 


20x50 


Cut Plane: 


45x50 


Outer Upper Shell: 


20x45 


Outer Lower Shell: 


20x45 


Upper Back Plane: 


45x50 


Lower Back Plane: 


45x50 



Of those domains, the sting, the wake, and the plane that was created to divide the 
upper and lower surface grids were included, identically, in both the upper and lower 
blocks to be joined as identical elements. It should be noted that, in this application, the 
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sting is the representation of the line forward of the waverider nose, created as an entity 
of the computational space. It does not refer to the sting used in wind tunnel tests to 
mount the vehicle in the tunnel test section. From these domains, along with the upper 
and lower surfaces of the body, the six-domain upper and lower 3D volumetric grids, 
known as blocks, were created. Inspection of each of the blocks using GRDDGEN3D, 
was used to determine the presence of negative cell volumes. Negative cell volumes 
resulted when the grid lines within the grid volume had tangled and the indices of the 
grid lines in certain areas were no longer sequential. The implication of this is that Ax or 
Ay used by the finite differencing scheme in OVERFLOW 1.8b will be a negative 
number, which is unacceptable to the code. 

The GRIDGEN3D elliptic solver, with a LaPlacian foreground control function, 
was used to resolve the negative cell volumes in both the upper and lower blocks. The 
LaPlacian function tended to drive the cells within the 3D block to equal volume cells 
throughout the total computational volume, thus those cells with negative volumes were 
driven to match the average cell volume, which was positive. While equal volume was 
unlikely to be achieved in a 3D grid containing a sharp body, applying this control 
eliminated the crossing of gridlines. It was possible to allow the elliptic solver to run for 
too many iterations. If this was done the number of negative cell volumes tended to 
increase. Periodically the elliptic solver was stopped, and the grid was examined to 
ensure the number of cells with negative volume were continuing to decrease. 

Regardless of the care taken not to over manipulate the grid, the elliptic solver tended to 
increase the distance between the body and its adjacent grid points at all points except the 
plane of symmetry, where grid points were held as fixed entities. 

FAST images of the upper and lower 3D blocks, following the elimination of the 
negative cell volumes, were examined to identify coordinate discrepancies. It was 
necessary to reverse the direction of the J-planes (the surfaces radiating from the body to 
the outer shell in lengthwise planes) and K-planes (surface planes stepping out from the 
body surface to the outer shell) of the upper block. This was accomplished using the 
GRIDED option to reverse plane indices. Following coordination of the plane indices of 



11 



the upper and lower blocks, the two blocks were joined, again using GRIDED. 
Additionally, one J-plane was added to the upper and lower back planes of the block to 
provide the planes of symmetry that facilitated the use of a central-differencing scheme in 
the computational phase of this analysis. This single block was then saved and examined, 
once again, for negative cell volumes using GRIDGEN3D. 

The figure below provides a perspective view of the 3-D computational block. At 
the center of this block is the surface of the waverider, seen as the black body. Final 
dimensions of the computational block were 41x45x50 grid points, for a total of 94,250 
points. 




Figure 3. Three-Dimensional Computational Block 
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Coordinate directions used by OVERFLOW1.8b were specified to have the J- 
direction increasing along the length of the body from nose to tail. The K-direction was 
specified as positive around the body, from the upper centerline, out to the wing tip, and 
back into the centerline of the lower surface. The L-direction is positive in the direction 
radiating out from the body to the shell of the 3D grid. 

B. SELECTED FLIGHT CONDITIONS 

To support the study of Reynolds number effects for both viscous and inviscid 
flow conditions, a range (0.3 < Moo < 6.0) of computations were conducted. A summary 
of the flight conditions for which analysis was performed is provided in Table 1. 



Mach Number Range 


Angle of Attack 
Sweep 


Engine Flow- 
Through 


0° AOA 


M« = 0.3 


0° AOA 


Inviscid 


Viscous 


Viscous 


Inviscid 


0.3 


0.3 


(degrees AOA) 


M« = 6.0 


0.4 


0.4 




0.5 


0.5 


-10 


0.6 


0.6 


-8 


0.7 


0.7 


-4 


0.8 


0.8 


-2 


0.9 


0.9 


0 


1.2 


1.2 


6 


2.0 


2.0 


2 


4.0 


4.0 


8 


6.0 


6.0 


10 



Table 1. Summary of Flight Conditions Analyzed 

Both viscous and inviscid flow solutions for the Price waverider were run for Moo 
= 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.2, 2.0, 4.0, and 6.0. Alpha sweeps were conducted at , 
Moo = 0.3. Additionally, the M» = 6.0 input file was modified to include boundary 
conditions simulating engine intake and exhaust features to examine the effect of engine 
flow-through at design conditions. OVERFLOWl .8b is optimized for M® = 5.0 flow 
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solutions. While the code performed well over a broad range of Mach numbers, the time 
required to compute low speed solutions and the decreasing stability of the code at Mach 
numbers greater than Moo = 5.0 began to limit its application. 

Figure 4 shows the number of iterations required to achieve convergence in the 
flow solutions at each Mach number. At the lower Mach number, and indeed even at M»o 
= 0.3, the computational time required to fully converge a solution is over 130,000 
iterations. This equated to over 50 hours of computational time, using all four processors 
of the Cray supercomputer. The Mach 0.5 solution required nearly 160,000 iterations due 
to the smaller time step that was required to keep the solution stable. The reason for the 
instability seen at M<» = 0.5 was not clear. At M«, = 6.0, the time step required to prevent 
the solution from becoming unstable was very small and the computational solution was 
particularly intolerant of even moderately aggressive attempts to speed up the 
convergence. 
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Figure 4. Number of Iterations vs. Mach Number 
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c. 



FLOW SOLUTION BATCH-FILE INPUT 



OVERFLOW!. 8b is a batch-run FORTRAN-based flow solver, which required 
the creation of input files to start each flow solution. Input files for each flight regime 
were created, using Reference 8, to reflect the flow conditions, desired turbulence model, 
computational controls and boundary conditions. Each of the batch input files was set up 
to apply a central-differencing scheme on a single grid computational space. Mach 
number and Reynolds number inputs reflected the desired freestream Mach number and 
the Reynolds number computed for the 15-inch model at sea-level conditions. Boundary 
conditions at the upper and lower back planes of the 3D grid were defined as planes of 
symmetry in the Y-direction. The wake cut was intended to be invisible to the flow and 
was defined as a fold-over cut flow-through plane. The outer shell was defined as a 
supersonic and/or subsonic inflow/outflow surface. 

Throughout the course of computations the input files were updated to reflect 
more aggressive time steps, allowable after a solution began to converge. The time steps 
were adjusted without regard to real-time. Sample input batch-files are provided in 
Appendix B. 

A summary of the general approaches used is given below. 

1. Subsonic versus Supersonic Flowfields 

In general, the allowable time steps for the finite difference calculations were a 
function of the freestream Mach number and maturity (i. e., how far the solution had 
converged) of the solution. Once the solution began to converge time steps could be 
increased, usually by factors of two until convergence. Allowable time steps were 
greatly reduced at M„ = 6.0 because of the velocity of the flow. The time step here had 
to be reduced by an order of magnitude, to dt = 0.001, to meet stability requirements. 
Tannehill, Anderson and Pletcher provide an excellent discussion of stability 
requirements and Courant-Friedrich-Levy (CFL) numbers [Ref. 12] . 

Transonic and supersonic flow solutions required the use of second-order 
smoothing in all three directions to accommodate discontinuities in the flowfield across 
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shocks. Shock discontinuities became present in the flowfield at freestream Mach 
numbers of ML» = 0.7 and higher. Prior to that Mach number, no special considerations in 
the input files were required. 

2. Inviscid versus Viscous Flowfields 

The input files for Euler (inviscid) solutions required fewer arguments as the flow 
description defined no viscosity effects in any direction nor was any turbulence model 
specified. Conversely the input batch-files for the viscous cases specified viscosity 
effects in the J- (around the waverider body) and L-directions (from the waverider body 
to the outer shell). Viscosity in the K direction (along the length of the waverider body) 
was neglected in order to reduce the computational time. Neglecting these streamwise 
viscous terms resulted in the solution of a parabolized form of the Navier-Stokes 
equations, which was computationally more efficient. The Baldwin-Lomax turbulence 
model was selected for use in all viscous cases. 

3. Engine Flow-Through Modeling 

Input batch-files for the engine flow-through for required additional boundary 
conditions that described the intake area as a region of extrapolated outflow (i. e., the 
flow approaching the intake was allowed to exit the computational domain through that 
region, shown in Figure 5). This figure shows the flow approaching the engine inlet from 
the right and penetrating the area of the inlet, seen as arrows inside the engine duct. 
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Figure 5. Engine Flow-Through Model Intake 



The engine exhaust are was defined as nozzle inflow (i. e., the mass flow, 
extrapolated from po and To, was injected into the computational domain at the exhaust 
plane). Engine exhaust flow was defined as freestream flow, shown in Figure 6. Here 
we see both the results of the external flow solution as green arrows and the 
superimposed engine flow-through exhaust, seen as the arrows in the freestream 
direction. 
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IV. RESULTS AND DISCUSSION 



A. SUBSONIC RESULTS 

Integral to the validation of computational solutions is the attainment of “steady- 
state” or convergence. As the boundary conditions propagate through the computational 
space, values at a given point are compared to the values of the previous iteration until 
the residual difference becomes negligible. In this study, a drop in residual values of 
three orders of magnitude was considered a converged solution. A sample convergence 
history is provided in Figure 7. 

To validate each of the flow solutions reported in this study, each flight condition 
solution was examined for convergence of the solution’s residuals, convergence of the lift 
and drag coefficients to a steady state, and representative surface pressure distribution 
and Mach number distribution within the flowfield surrounding the vehicle. 

1. Mach 0.3 Inviscid Results 

The significant computational time required to conduct low subsonic 
computational analysis limited low-subsonic solutions to a minimum Mach number of 
Mac = 0.3. Wind tunnel tests provided data up to a maximum Mach number of M* = 

0.15. While the computational analysis provided here is not entirely coincidental with 
existing wind tunnel data, the solution obtained demonstrated fair representation of the 
low speed flow characteristics of the Price waverider. 

Examining the M* = 0.3 inviscid solution, shown in Figure 7, shows that the 
convergence history achieved just over three orders of magnitude convergence in 126,200 
iterations. The short period oscillations were the result of periodic vortex formation and 
subsequent shedding at the trailing edge of the vehicle. Also present in the convergence 
history was a long period oscillation occurring approximately once every thousand 
iterations. This was speculated to be the result of acoustic waves reverberating within the 
computational space. While no direct evidence within the solution of the flowfield 
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Residual 



supported that conclusion, their presence in the subsonic solutions only supported the 
conclusions that an acoustic phenomena of some type had occurred. 




Figure 7. Mach 0.3 (Inviscid) Residual Convergence History 



Figure 8 shows the drag coefficient history. Note that the drag coefficient also 
converged to a steady state as the solution converged. The final achieved drag coefficient 
was computed was Cd = 0.0194. 
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Figure 8. Drag Coefficient vs. Number of Iterations for M» = 0.3, Inviscid 

Figure 9 shows the lift coefficient history. This force coefficient converged to a 
stable value coincident with the drop in residuals. At M* = 0.3, the lift coefficient 
converged to Cl = -0.03 10. The presence of a negative lift coefficient at subsonic 
speeds was the result of the body geometry that produced lift at supersonic speeds. While 
the wedge angle of the vehicle forebody created a mean chord line, that was a straight line 
at a positive angle of attack, the nozzle expansion aft of the engine cowl and the shape of 
the lower wing surface created a negative camber in the vehicle, giving rise to negative 
lift at zero degrees angle of attack. 



21 




Figure 9. Lift Coefficient vs. Number of Iterations for M*. = 0.3, Inviscid 



Using FAST images, the flowfield solution was graphically examined. Figure 10 
shows the distribution of the normalized pressure (pressure is normalized to the 
ffeestream static pressure) on the body surface and the Mach number contours normal to 
the body in the surrounding flowfield. The normalized pressure was equal to the 
ffeestream pressure over the entire upper surface. In the inviscid case, the upper surface 
flow was unperturbed, as the waverider upper surface was aligned with the ffeestream. 
The lower surface pressure distribution illustrated a slightly higher pressure area on the 
lower surface of the vehicle forebody than on the lower aftbody, indicating a potential for 
significant pitch-up in the vehicle at low speeds. This was borne out by the computed 
pitching moment coefficient of Cm = -0.033. Neither region displayed a pressure 
distribution that was capable of providing lift, as was confirmed by the negative lift 
coefficient. 
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Figure 10. Normalized Surface Pressure and Mach Number Contours, Moo = 0.3, Inviscid 



The Mach number contours confirmed the nearly unperturbed upper surface, 
which contrasted with the large Mach number gradients at the trailing edge and the less 
extreme but critical, Mach number contours along the lower surface. The Mach number 
contours at the nose of the vehicle indicated that the stagnation point was located on the 
lower surface of the body slightly aft of the nose of the vehicle. Contours along the lower 
surface outlined the flow compression at the beginning of the first and second engine 
inlet ramps and subsequent expansion as the flow turned back toward the ffeestream at 
the engine cowl and upward again at the engine nozzle. Further away from the body, grid 
spacing was less dense, which resulted in the Mach contour lines not being represented as 
smooth curves. This lack of smooth contour lines, indicated a grid dependency within the 
solution. 
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2 . 



Mach 0.3 Viscous Results 



Viscous terms were added to the partially converged inviscid solution at 66,200 
iterations in order to capitalize on the convergence already achieved during the inviscid 
computation. OVERFLOW 1.8b allowed the viscous modification to previously inviscid 
solutions. This addition had several notable impacts. Shown in Figure 11, residuals rose 
sharply as the flow was “virtually” perturbed by the addition of viscosity, spiking at well 
over the initial maximum, but recovering within 1000 iterations. As is seen in Figure 11, 
short period oscillations remained present, as the more physically representative wake in 
the viscous model continued to generate and shed vortices. The sonic wave reverberation 
persisted as well. It should be noted that the addition of viscous terms enabled a slightly 
more aggressive time step to be used in the remaining computations. 

It was hoped that the addition of these viscous terms would serve to stabilize the 
solution and gain faster convergence, while it can be seen that the opposite result was 
achieved. Despite the more aggressive time step and expectations of a more rapid 
convergence, the slope of the declining residuals when compared with the inviscid case 
appeared to be the same. 
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Figure 1 1. Mach 0.3 (Viscous) Residual Convergence History 



Figures 12 and 13 illustrate the impact of the introduction of viscous terms on 
both the lift coefficient and the drag coefficient. The computation of a more physically 
realistic wake, enhanced the vehicle performance in both axes. A small, but notable, drop 
in drag occurred with the addition of viscosity, which demonstrated that the addition of 
viscous drag over the body was more than compensated by the significant reduction in 
pressure drag achieved as the viscous terms forestalled the onset of separation at the 
trailing edge. Similarly, the drag coefficient, which continued to oscillate, shifted 
upwards to a new mean lift coefficient at the introduction of the viscous terms. 
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Figure 12. Drag Coefficient vs. Number of Iterations for M* = 0.3, Viscous 
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The graphical representation of the surface pressure and Mach number contours in 
the viscous case at — 0.3 are shown in Figure 14. The presence of Mach number 
variation in the flowfield over the upper surface was an indication of the boundary layer 
development and propagation into the flowfield. The lower surface boundary layer was 
visible along the engine cowl. The region of the trailing edge showed a narrower wake 
that was believed to be more representative of the expected physical wake, due to the 
presence of viscosity in the flow. The computed viscous wake (Figure 14) was smaller 
and better defined than the wake in the inviscid case (Figure i 0), which indicated less 
separated flow in the region of the trailing. 
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3. 



Subsonic Trends 



The convergence histories, force coefficient histories and FAST images for all 
solutions can be seen in Appendix C. Flow solutions for the remainder of the subsonic 
flight conditions demonstrated many similar characteristics. Residual plots demonstrated 
the requisite convergence to a steady state solution. The addition of viscous forces after a 
number of inviscid iterations caused a spike to appear in the residual files similar to the 
spike seen in Figure 11. Beginning at M® = 0.5, the presence of the long-period acoustic 
resonance dampens out in the early iterations leaving only the shot period vortex- 
shedding oscillation evident within the residual plot. An additional undamped oscillation 
occurred in the nearly converged viscous solutions at M® = 0.7 (See Figure C-10) , 0.8 
and 0.9. Visual inspection of the solutions of these flight conditions using FAST showed 
the separation and recirculation of flow in the wake of the vehicle. This separation was 
believed to be the cause of the unsteadiness in the residual. 

With increasing Mach number the stagnation point, evidenced by the Mach 
contours seen at the nose, migrated forward from a position aft of the nose to a position 
coincidence with the nose of the waverider (See Figure C-47 to Figure C-60). 

Lift and drag coefficient histories showed very little difference in patterns across 
the subsonic Mach number spectrum; each file showed oscillating values which damped 
to a single value as the solution converged. Behaviors of the lift and drag coefficients, 
with the addition of viscous terms, were similar to the case of M® = 0.3 viscous flow. 
Trends in the values of these coefficients are the subject of the later section on Reynolds 
number effect. 

B. SUPERSONIC RESULTS 

The supersonic solutions had the advantage of requiring significantly less 
computational time to obtain a fully converged solution. The Moo = 4.0 solution 
converged in 4,100 iterations, as opposed to the 160,000 iterations required for Moo = 0.5. 
In addition to the speed at which the codes ran, a number of other differences between the 
subsonic and supersonic solutions exists. 
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1 . 



Mach 6.0 Inviscid Results 



Figure 13 shows the convergence history for the inviscid Moo = 6.0 solution. This 
solution required some 13,000 iterations due to the small time step required to prevent the 
solution from going unstable. Nonetheless, the computational time was an order of 
magnitude lower than the subsonic solutions. Note that both the long-period oscillations, 
attributed to acoustic wave reverberation within the computational space, and the short 
period oscillations, attributed to trailing edge vortex shedding, seen in Figures 7 and 11, 
were absent in the M® = 6.0 convergence history. 

The shape of the convergence history, when compared to the subsonic results 
(Figures 7 and 1 1) , was markedly different. The residuals failed to drop off at the nearly 
steady rate seen in the subsonic convergence histories, but, instead, maintained a nearly 
constant value until late in the computation and finally dropped off sharply. The reason 
for the convergence history remaining nearly constant and then suddenly dropping off 
was that, within the flowfield of the waverider at supersonic conditions, the bulk of the 
flowfield was largely unaffected, thus the residuals did not converge until the effected 
sections were fully updated. Once the solution began to update the region close to the 
body, the solution converged rapidly. 
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Figure 15. Mach 6.0 (Inviscid) Residual Convergence History 



Despite the lack of a fully converged solution early in the computation the drag 
(Figure 16) and lift (Figure 17) coefficients rose quickly and remained nearly steady 
throughout the computation. Also noteworthy was the drop in the drag coefficient by 
half, to Cd = 0.00439 from the = 0.3 inviscid solution of Cd = 0.0194. As expected, 
the lift coefficient at M« = 6.0 is positive. At Mach 6.0, the converged lift coefficient 
was Cl = 0.018. 
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Figure 17. Lift Coefficient vs. Number of Iterations for = 6.0, Inviscid 
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Figure 18 shows the normalized pressure on the surface of the waverider as well 
as the computed Mach number contours in the flowfield surrounding the body. The 
pressure differential between the fore and aft bodies of the vehicle appeared less 
pronounced than at = 0.3, however the scale of the legend was changed to 
accommodate the large pressure rise at the engine intake. Pressures of nearly 14 times 
the freestream static pressures were seen in this region of the flow. 




Figure 18. Normalized Surface Pressure and Mach Number Contours, Moo = 6.0, Inviscid 

The Mach contours illustrate the formation of shocks at the nose of the vehicle 
and at the engine inlet. From the NACA Report 1135 Oblique Shock Chart [Ref. 13], the 
predicted shock angle for M® = 6.0 flow over the 5.2° wedge of the vehicle nose is 13°. 
Although difficult to measure from the FAST image, this appeared to be the angle at 
which the shock was computed. At that freestream Mach number the shocks over the 
engine inlet two-ramp configuration coalesced into a single oblique shock. The engine 
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inlet however, when modeled as a solid body, dramatically changed the flow angle and 
caused a distinct shock to form at the engine inlet. Here, where the engine inlet was 
modeled as a solid body, the engine inlet shock formed at a much higher angle than the 
shocks off the nose of the vehicle. 

2. Mach 6.0 Viscous Results 

Due to the tendency of the = 6.0 solution to go unstable with minor changes 

in the inputs, the viscous solution was computed from zero iterations. This convergence 
history (Figure 19) shows that the viscous solution converged more slowly than the 
inviscid M® = 6.0 solution (Figure 15). The viscous effects required additional iterations 
to allow development of the boundary layer. The drag coefficient and lift coefficient 
histories shown in Figures 20 and 21, also reflected the presence of viscous effects from 
the beginning of the computations, rather than showing a change within their values when 
viscous forces were added. 




Figure 19. Mach 6.0 (Viscous) Residual Convergence History 
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Figure 20. Drag Coefficient vs. Number of Iterations for M» = 6.0, Viscous 




Figure 21. Lift Coefficient vs. Number of Iterations for M* = 6.0, Viscous 
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In Figure 22, the presence of Mach contours in the flowfield near the upper 
surface of the Price waverider was evidence of the boundary layer formation on that 
surface. Shocks were more smeared than in the inviscid case and the wake region 
extended over a larger distance downstream of the trailing edge. 




3. Supersonic Trends 

Inviscid and viscous solutions were computed for M«> = 1.2 and M«. = 4.0, with 
similar results to the M» = 6.0 solutions. The engine flow-through model was only 
computed for the M*> = 6.0 flight condition. The results for additional supersonic cases 
are provided in Appendix C. 
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c. 



ENGINE FLOW-THROUGH MODELING 



The presence of the extremely high pressure forward of the intake, seen in Figures 
18 and 22, did not necessarily represent the true nature of the flowfield of an operating 
vehicle. To provide data on a more representative configuration engine flow-through was 
modeled as a follow-on test case to the = 6.0 computations. This configuration was 
not modeled for the subsonic case, although the impact is expected to be less significant 
because the flow in front of an engine in subsonic flight is expected to be slowed 
significantly at the diffuser and thus will present a surface that more closely resembles a 
solid body. 

The Mach 6.0 engine flow-through computation could not be run from the 
previously converged = 6.0 solution because it not only required modification to the 
boundary conditions, but required the definition of a new set of boundaries. Flow- 
through was extrapolated from the pressure forward of the intake boundary and the 
exhaust was modeled as freestream pressure extrapolated flow. 

Convergence history for the M„ = 6.0 engine flow-through case (Figure 23), 
shows more moderate drops in the residuals than in the initial = 6.0 solution (Figure 
15). Although the residuals began to drop off between 7000 and 8000 iterations, as seen 
in the previous Mach 6.0 solution, the engine Flow-Through model required nearly twice 
as many iterations to converge, as the pressure extrapolations from previous iterations 
were computed as input to follow-on iterations. 
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Figure 23. Mach 6.0 (Inviscid) Engine Flow-Through Convergence History 

The resulting solution was computed to show decreases in both lift and drag 
coefficients over the M* = 6.0 solution. The lift coefficient, calculated in the inviscid 
M® = 6.0 solution was Cl = 0.018, was found in the engine flow-through model to be 
Cl = 0.016. The drag coefficient was found to be Cd = 0.003 for the engine flow-through 
model, whereas previous M» = 6.0 solution (without engine flow-through) computed the 
drag coefficient to be Cd = 0.0044. In this calculation, the declining lift was outpaced by 
the improvements in the drag, resulting in better overall L/D performance than the 
original calculations demonstrated. 

As is illustrated in Figure 24, the pressure at the intake dropped to only twice the 
freestream static pressure, as opposed to the factor of 14 seen in the initial MU = 6.0 
solution. Mach number contours also showed that the resulting shocks occurred closer to 
the body. 
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D. REYNOLDS NUMBER EFFECTS 

Performance of the Price waverider across the spectrum of flight Mach numbers 
was compared using the lift and drag coefficients and the overall performance factor of 
the lift-to-drag ratio. The values for lift and drag coefficients were tabulated from the 
OVERFLOW1.8b output file containing force and moment data. The force and moment 
file provided the force and moment contributions due to pressure and friction. These 
results are compiled in Tables 2 and 3. In the inviscid case the contributions due to 
friction were zero and have not been reported in Table 2. In the viscous cases, the 
computed contribution of the viscous effects to the lift and drag coefficients, were, in 
nearly all cases, rounded to zero. 
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Mach 

Number 


C L 


Cd 


L/D 


Cm 


0.3 


-0.0310 


0.0194 


-1.60 


-0.0387 


0.4 


-0.0329 


0.0171 


-1.92 


-0.0405 


0.5 


-0.0344 


0.0160 


-2.15 


-0.0424 


0.6 


-0.0377 


0.0159 


-2.37 


-0.0460 


0.7 


-0.0402 


0.0162 


-2.48 


-0.0494 


0.8 


-0.0421 


0.0178 


-2.37 


-0.0530 


0.9 


-0.0308 


0.0230 


-1.34 


-0.0463 


1.2 


0.0078 


0.0302 


0.26 


-0.0167 


2.0 


0.0312 


0.0152 


2.06 


0.0118 


4.0 


0.0224 


0.0064 


3.50 


0.0108 


6.0 


0.0184 


0.0044 


4.20 


0.0094 



Table 2. Summary of Inviscid Solution Force Coefficients 



Mach 

Number 


Cl 

pressure 


C L 

friction 


Cl total 


Cd 

pressure 


Cd 

friction 


Cd total 


L/D 


Cm 


0.3 


0.0239 


0.0000 


-0.0238 


0.0169 


0.0001 




-1.40 


-0.0332 


0.4 


0.0203 


0.0000 


H&O203 


0.0156 


0.0001 


..... . 7 . 


-1.29 


-0.0306 


0.5 


0.0176 


0.0000 


m$m 


0.0150 


0.0001 




-1.17 


-0.0289 


0.6 


0.0149 


0.0000 


-0.0149 


0.0149 


0.0001 


13M 


-0.99 


-0.0272 


0.7 


0.0109 


0.0000 


-0.0109 


0.0152 


0.0001 




-0.71 


-0.0247 


0.8 
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Table 3. Summary of Viscous Solution Force Coefficients 



These results are compiled for both the viscous and inviscid solutions and shown 
graphically. Figure 25 shows the lift coefficient as a function of Mach number. The 
inviscid subsonic lift behavior showed the detrimental effects of separation, although it is 
not clearly understood what specific effect caused the inviscid lift coefficients to diverge 
so dramatically from the viscous lift coefficients subsonically. In the supersonic region. 
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viscous and inviscid lift coefficients showed very good agreement with each other. The 
implication here is that inviscid solutions closely approximate the behavior of viscous 
flows, and as such, the viscous effects may be negligible in the initial design phase of 
supersonic vehicles. The viscous solutions consistently computed higher values for the 
lift coefficient, which demonstrated the beneficial effects of the attached flow. A 
significant feature of this configuration was that the lift coefficient does not become 
positive until the flow is supersonic. In short, the Price waverider, in its current 
configuration, produced no positive aerodynamic lift until the freestream Mach number 
became sonic. Studies of the drag coefficient and pitching moment are required to 
determine the impact of not have positive lift at zero degrees angle-of-attack. 

Wind tunnel testing initially performed by Cedrun [Ref. 6] provided a basis of 
comparison for Huffs follow-on testing [Ref. 7]. Huff demonstrated that his data 
duplicated Cedrun’s original results. Due to the formatting of Cedrun’s work, the lift and 
drag coefficients, with wall and blockage corrections for the wind tunnel, were not 
provided in tabular form, so for greater accuracy, the values from Huffs work are shown 
in the comparisons that follow. Since the data of both sets of testing agree, the 
comparisons made here to Huffs data can be applied to Cedrun’s work as well. 

The lift coefficient from Huffs highest speed wind tunnel tests is plotted on 
Figure 25 at = 0. 15. Seen here it lies slightly below the extrapolated path of the 
viscous solution of computed lift coefficients, demonstrating that the wind tunnel test 
results are slightly less optimistic than the computed solution. The reason for this is 
believed to be the addition of the trailing edge cusp on the computational waverider 
model. 

The single engine flow-through configuration lift coefficient is also plotted on 
Figure 25 and shows the slight decrease in lift at = 6.0. This is due to the angle of the 
shock formed at the engine inlet. In the engine flow-through model, this shock was 
slightly more oblique than the shock seen in the solid body engine model at = 6.0. 

The more oblique shock in the engine flow-through model causes a smaller pressure rise 
across the shock, thus the pressure on the lower surface of the waverider body is less. 
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In Figure 26, the plot of drag coefficient versus Mach number, we see the 
transonic drag rise across = 1.0 in both the viscous and inviscid solutions. The 
dominant effect of the pressure drag over any viscous forces is supported by the close 
agreement of these two sets of data through the entire Mach number range computed, 

0.3 < M„ < 6.0. The trough in drag in the subsonic region was an unexpected occurrence. 
The subsonic flowfield results, shown in Appendix C, indicate that the stagnation point 
migrated from a position slightly aft and under the nose of the vehicle at M*. = 0.3, to a 
position at the nose for M» = 0.5 or 0.6. This stagnation point migration presented a 
more streamlined body to the freestream flow and resulted in a decrease in the drag 
coefficient. This streamlining effect was overcome by the drag rise in the transonic 
region. This transonic drag rise was the result of the presence of regions of supersonic 
flow that had to shock back to freestream along the surface of the body. The presence of 
these shocks caused the early onset of separation, which increased the pressure drag felt 
by the body. 
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Figure 25. Lift Coefficient vs. Mach Number 
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Included within the graphical representation of drag coefficients is the M* = 0.15 
value of Cd from Huffs wind tunnel tests. This point lies slightly above the 
extrapolation of the computed solutions. Also plotted is the drag coefficient for the M» = 
6.0 engine flow-through configuration. As in Figure 24, the modeling of engine flow- 
through reduced the shock angle, and presented a more streamlined body to the 
freestream flow, which reduced the drag significantly. 




Mach Number 

Figure 26. Drag Coefficient vs. Mach Number 

Overall performance of the Price waverider is shown in Figure 27. The plot of 
lift-to-drag ratio versus Mach Number shows the optimum performance at M® = 6.0. 
This result was expected, since the vehicle configuration was optimized for M«, = 6.0, 
although this performance indicator did not reflect the product of lift-to-drag ratio and 
specific impulse. This study did not include the estimation of engine performance. It 
should also be noted that, for purposes of consistency throughout the study, all Mach 
numbers and Reynolds numbers were computed at sea level, an unrealistic environment 
for the M» = 6.0 flight condition. 
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The engine flow-through performance indicated that the decrease in drag in this 
configuration outpaced the decline in lift performance. Also plotted are the Moo = 0.15 
lift-to-drag ratio for Huffs wind tunnel test data which lies in the plane of the 
extrapolated L/D curve from the computed solutions. 



12 

10 

5 8 

cz 

* 6 
DJD 
zz 

5 4 

2 2 

£ 

3 0 
-2 
-4 









i i i ; 

♦ ■ L/D - Inviscid 

* L/D - Viscous 

L/D - with Engine Through flow 
— — L/D - Kuchemann’s Design Target 
— X— Wind Tunnel Data (Huff) 






\ 
































~ ' ” - — 









^=4 














m 





































0 1 2 3 4 5 6 

Mach Number 



Figure 27. Lift-to-Drag Ratio vs. Mach Number 



Maximum L/D performance for a given design Mach number waverider is shown 
for the supersonic range. Shown here, Kuchemann’s maximum on-design performance is 
compared to the off-design performance data for the Price waverider. This estimate was 
based on Kuchemann’s equation [Ref. 3] of what should be attainable in optimized 
design: 

I _ 3(A/,+3) 

D M„ 

Kuchemann’s equation is based on empirical data and predicts that, in the limit as 
Moo goes to infinity, L/D will asymptote to 3.0. At the design Mach number of M» = 6.0, 
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all L/D performance values are shown to converge, demonstrating close agreement 
between the computed and the theoretical L/D performance of the vehicle. The predicted 
value of L/D at M«>= 6.0 is 4.5. This performance is closely matched by the computed 
values of L/D = 4.19 for the inviscid solution, L/D = 4.54 for the viscous solution and 
L/D = 5. 1 7 for the engine flow-through solution. Note that both the viscous solution and 
the engine flow-through solution exceed that maximum predicted value for L/D. 
Anderson [Ref. 1], disputed Kuchemann’s performance limits, estimating that L/Dmax 
for supersonic vehicles should be about 50% higher than Kuchemann’s predictions for 
Mach 6.0 optimized vehicles. Anderson, however, based that dispute on a class of 
waverider-configured vehicles generated using his viscous-optimization routine. While 
the viscous effects have been taken into account in the computations reported herein, the 
Price waverider was designed using an inviscid optimization routine. 

Figure 28 shows the value of the pitching moment coefficient for the computed 
solutions. Based on the coordinate system used, the positive pitching moment represents 
a nose down pitching force. Without employment of control surfaces, the Price 
waverider appears to be statically unstable. The datum employed here is the y-axis, the 
nose of the vehicle. The inviscid solution shown here continues to demonstrate the 
adverse effect of the non-physically representative wake in the CFD solution, and the 
subsonic values for Cm do not reflect a realistic solution. 
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Figure 28. Pitching Moment Coefficient vs. Mach Number 



As is seen in Figures 25 through 28, the computed force coefficients demonstrated 
agreement with the experimental data from Huffs low-speed wind tunnel testing. 

Overall, the modeling of engine flow-through improved the performance data. Close 
agreement between the viscous and inviscid solutions in the supersonic region supported 
the premise that, for a first approximation, viscous effects can be neglected in the initial 
design phases. This argument cannot be supported for the subsonic region, as the effects 
of the non-physically representative wake cause lift, LTD performance and pitching 
moment parameters to diverge from the viscous case. 

E. ALPHA SWEEP 

Computations for the range of angles of attack of -10° < a < +10° were 
conducted at M«, = 0.3, for viscous conditions. The convergence history, as well as lift 
and drag coefficient history, were similar to those seen in the zero degree angle of attack 
solution for = 0.3, viscous flow. The calculated effect on the lift coefficient is shown 
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in Figure 29. Wind tunnel data for M® = 0. 1 5 from HufFs thesis are shown for 
comparison. The computed values demonstrated a lift curve slope of Cl g = 0.884, 
compared to the wind tunnel test lift curve slope of Cl« ~ 2.526. This discrepancy is 
most likely due to the coarseness of the grid. It is also possible that a better match may 
be obtained by using a different turbulence model for the computations. Cedrun’s test 
results, although not plotted, matched Huffs data. Cedrun noted no marked difference in 
the lift curve slope relative to Mach number, so the impact of the difference in freestream 
Mach number between the wind tunnel data and the computed solutions was discounted 
as the cause for the differential in lift curve slope. 




Figure 29. Lift Coefficient vs. Angle of Attack 

The onset of stall was not predicted to occur within the range of angles of attack 
for which solutions were computed. Cedrun’s lowest speed wind tunnel test (M® * 
0.003) concluded that stall occurred at 34 degrees AO A. Neither Huff nor Cedrun 
achieved stall in higher speed wind tunnel tests, due to the total force vector exceeding 
the wind tunnel sting force constraint. 
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Figure 30 shows the impact of angle of attack on the drag coefficient for both 
wind tunnel data and computed solutions. The computed drag coefficient was 
consistently lower than the wind tunnel test data The optimistic prediction of drag 
coefficients was presumed to be the result of the addition of the trailing edge cusp. The 
addition of the cusp also caused a smaller variation in drag coefficient over the range of 
angles of attack, thus flattening the curve when compared to the wind tunnel data. 
Minimum drag was found to occur at zero degrees angle of attack. 




-10 -5 0 5 10 

Angle of Attack 

Figure 30. Drag Coefficient vs. Angle of Attack 

A final comparison was made between the wind tunnel data and the computed 
solutions' L/D performance. Shown in Figure 3 1, the measurements on the wind tunnel 
model showed a steeper variation of L/D with changing angle of attack. The lift curve 
slope variation was the largest contributor to the difference in L/D performance, however 
the addition of the trailing edge cusp also played a role in the differences in performance. 
Maximum L/D performance was calculated to be L/D = 2.87, which occurred at eight 
degrees angle of attack. The wind tunnel data showed a maximum of L/D = 3.58, which 
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occurred at four degrees of angle of attack. This value may have continued to rise had 
higher AOA’s been tested. At negative angles of attack the computed L/D performance 
predicts a minimum L/D between -4 and -8 degrees, closely matching the wind tunnel 
measurements, which measured a minimum L/D at -6 degrees. 




Angle of Attack 

Figure 3 1 . Lift-to-Drag Ratio vs. Angle of Attack 

The normalized surface pressure and Mach number contours for the positive 10 
degrees angle of attack case are shown in Figure 32. The sharp Mach number contour 
gradient showed the presence of the boundary layer on the upper surface of the body. 
Mach contours outlined the stagnation point located aft of the nose on the lower surface 
and indicated a region of separated flow on the upper surface. The portion of the upper 
surface where flow was separated was also indicated by the region of lower surface 
pressure seen along the leading edge from the nose to mid-chord. 
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F. CFD VALIDITY 

Based on the data presented above, the CFD appeared to calculate solutions, 
which represented realistic flowfields. While computational limitations prohibited 
calculation of the flowfield solution below M* = 0.3, the computed lift and drag 
coefficients and L/D for 0.3 < M«> < 6.0 trended toward the values obtained in the Moo = 
0.15 tests conducted by Huff and Cedrun. While precise comparison was not made, due 
to the mismatch of Mach numbers in the CFD and wind tunnel test, it appeared that the 
CFD trends were slightly optimistic. The optimism of the CFD calculations for lift and 
drag was believed to be the result of the addition of the trailing edge cusp to the CFD 
model. 

Angle of attack sweeps showed disagreement in all computed performance 
parameters. While some of the discrepancy was attributable to the addition of the trailing 
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edge cusp, the significant difference in the lift curve slope is believed to be the result of 
the grid coarseness and turbulence model. Despite disagreement at both positive and 
negative angles of attack, the solutions at zero degrees angle of attack were in close 
agreement. For that reason, the zero degree AOA data were still considered valid. 

The flight conditions considered herein represent the Mach number and Reynolds 
number calculated for the 15-inch model at sea level. To better approximate the 
boundary layer development, Reynolds numbers should be chosen to reflect the length of 
the Price waverider, rather than the model root mean chord. 



G. CFD LIMITATIONS 

For purposes of this study, a broad range of Mach numbers was examined, 
however the extent to which existing data could be matched was limited by the flow 
solver. Flow solutions were not obtained for Mach numbers below M» = 0.3, whereas 
wind tunnel test data was not available above M® = 0.15. This mismatch prevented the 
validation of CFD results with corresponding wind tunnel data. 

Difficulty in creating a grid that correctly modeled the blunt trailing edge, led to 
the addition of the trailing edge cusp to the CFD model. This difference in configuration 
is believed to have caused the lift and drag coefficients to be more optimistic than the 
wind tunnel test data. 

The use of the elliptic solver, GRIDGEN3D, to resolve the negative cell volumes 
in the grid also tended to pull grid points away from the body of the waverider. Due the 
migration of grid points away from the body of the waverider, the boundary layer 
calculations, while seen in the Mach number contours, may not be representative of 
actual boundary layer size. 
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V. 



CONCLUSIONS 



The computational fluid dynamic analysis of the Price waverider was 
demonstrated to produce valid results for subsonic and supersonic flowfield 
computations. The prediction of lift and drag coefficients and L/D closely resembled the 
data obtained in wind tunnel testing. Although the lift and drag coefficients appeared 
slightly more optimistic than the wind tunnels tests, this was likely the result of the 
addition of the cusped trailing edge. While that modification was made to facilitate the 
grid generation, the physical impact was to reduce the pressure drag and increase the lift 
through the expansion over the upper surface trailing edge. The differences in the 
computed and measured values were in the same order of magnitude as the tunnel wall 
corrections applied to the wind tunnel results. 

The use of computational methods provided the opportunity to expand the body of 
knowledge available defining the behavior of flow surrounding the Price waverider over 
a broad range of flight conditions. 

The angle of attack sweep calculation, while showing similar trends in 
performance parameters, did not achieve very good agreement, except at zero degrees 
angle of attack. The lift curve slope showed the largest discrepancy when compared to 
the wind tunnel data. This was believed to be the result of the grid coarseness. Further 
study is required to determine whether grid refinements or use of a different turbulence 
model would improve agreement. 

Instability of the solution at high Mach numbers precluded the extension of CFD 
analysis to flight conditions above the design condition. Time steps required to maintain 
a stable solution at M„ = 6.0 were an order of magnitude smaller than any other 
computation. Similarly, the excessive computational time required to compute the low 
subsonic solutions, prohibits the use of OVERFLOW1.8b to speeds lower than = 0.3. 

The computed flow solutions demonstrated several areas of concern. The first 
concern was the lack of positive lift generation in subsonic flight at zero degrees angle of 
attack. With a tendency to pitch up at this angle of attack present, the risk of flying at a 
slightly nose up attitude, to gain positive lift, is that this attitude will exacerbate the 
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pitching moment. The second area of concern was the transonic drag rise and the effect 
of the drag trough in the subsonic regime on the specific excess power curves. 

A number of opportunities exist to improve the body of knowledge regarding the 
Price waverider and in the extension of the existing analysis. 

A. GRID REFINEMENTS 

The existing grid provided a computational space that allowed the computation of 
reasonable solutions, however, there were a number of areas that needed improvement. 
The existing grid, after using the elliptical solver to resolve the presence of negative cell 
volumes, failed to maintain the initially specified grid spacing near the body except at the 
centerline where all points were held as fixed points. This spacing contributed to 
inaccuracies in the computation of viscous drag insofar as the velocity gradient in the 
boundary layer was dependent on the existence of sufficient points near the body to 
capture the gradient correctly. Additionally, the lift curve slope in the alpha sweep did 
not mirror the wind tunnel results. Grid density was believed to be the reason for this 
disagreement. 

The addition of the cusp at the trailing edge of the Price waverider, while a likely 
design change in the real world, did not reflect the vehicle as-tested. It was seen that the 
lift and drag forces computed over the existing grid were optimistic when compared to 
the wind tunnel test data. Use of multi-block grid generation techniques may improve the 
opportunity to match the grid geometry to the wind tunnel model 

It should be noted that the leading edge of the computational model was 
represented by a single point line, with no leading edge radius. This contributed to a 
departure from both the model, as-tested, and the real-world vehicle requirements. While 
the leading edge radius of the Price waverider was not included in the design criteria, it 
can be safely assumed that the leading edge would have to have some bluntness. An 
infinitely thin leading is unrealistic to build, and would not be able to withstand the high 
heating conditions produced at hypersonic speeds. This lack of bluntness at the leading 
edge contributed to unrealistic boundary layer development. 
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In general, the grid was fairly coarse, as only 45 points along the length of the 
waverider body were used during the computations. Grid dependence could be seen in 
the unevenness of the Mach contours in several of the solutions. 

B. OPPORTUNITIES FOR FURTHER TESTING 

The use of a multi-grid computational model would enable the closer examination 
of some aspects of the flowfield. This method could also be used to make modifications 
to the basic waverider, to include cockpit, vertical tails or high-lift devices to gain a better 
understanding of the aerodynamic impacts of varying configurations. Included in this list 
is the modeling of the engine as a thrust producing element rather than as a solid body or 
a flow-through duct. As with the employment of a cusped trailing edge, a number of 
configuration changes or proposed changes will require analysis before the feasibility of 
the Price waverider, or a derivative, can be fully assessed. 

This study did not examine the effects of yaw angles or the combinations of 
angle-of-attack and yaw on the flowfield. The opportunity, therefore, exists to further the 
understanding of the flight characteristics of the Price waverider at those off-design flight 
conditions. 

To date, supersonic wind tunnel testing has not been accomplished for the Price 
waverider. A study that combined supersonic wind tunnel test data, with the results 
provided herein would serve to validate the CFD code further. 
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APPENDIX A. WAVERIDER SURFACE GRID POINTS 
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0.76178 


0.36158 


0.06792 


0.76178 


0.36934 


0.05668 


0.76178 


0.37524 


0.06929 


0.76178 


0.38451 


0.06264 


0.76178 


0.38860 


0.07109 


0.76178 


0.38642 


0.06342 


0.76178 


0.39026 


0.07134 


0.76178 


0.39951 


0.06906 


0.76178 


0.40162 


0.07329 


0.76178 


0.41431 


0.07592 


0.76178 


0.41431 


0.07592 


34 






34 






0.81896 


0.00000 


0.00000 


0.81896 


0.00000 


0.10925 


0.81896 


0.01629 


0.00083 


0.81896 


0.01870 


0.10925 


0.81896 


0.03258 


0.00172 


0.81896 


0.03735 


0.10925 


0.81896 


0.04887 


0.00264 


0.81896 


0.05595 


0.10925 


0.81896 


0.06515 


0.00360 


0.81896 


0.07450 


0.10925 


0.81896 


0.06824 


0.00381 


0.81896 


0.07801 


0.10925 


0.81896 


0.06839 


0.00382 


0.81896 


0.07818 


0.10925 


0.81896 


0.06854 


0.00383 


0.81896 


0.07828 


0.07884 


0.81896 


0.08143 


0.00465 


0.81896 


0.08959 


0.07862 


0.81896 


0.09770 


0.00585 


0.81896 


0.10727 


0.07815 


0.81896 


0.11397 


0.00709 


0.81896 


0.12483 


0.07762 


0.81896 


0.13022 


0.00849 


0.81896 


0.14225 


0.07694 


0.81896 


0.14647 


0.00999 


0.81896 


0.15952 


0.07620 


0.81896 


0.16269 


0.01166 


0.81896 


0.17663 


0.07542 


0.81896 


0.17890 


0.01350 


0.81896 


0.19356 


0.07463 


0.81896 


0.19509 


0.01552 


0.81896 


0.21029 


0.07382 


0.81896 


0.21125 


0.01776 


0.81896 


0.22681 


0.07305 


0.81896 


0.22736 


0.02022 


0.81896 


0.24310 


0.07235 


0.81896 


0.24345 


0.02298 


0.81896 


0.25914 


0.07169 


0.81896 


0.25949 


0.02599 


0.81896 


0.27494 


0.07119 


0.81896 


0.27376 


0.02891 


0.81896 


0.28884 


0.07087 


0.81896 


0.27545 


0.02928 


0.81896 


0.29049 


0.07088 


0.81896 


0.29136 


0.03293 


0.81896 


0.30574 


0.07073 


0.81896 


0.30718 


0.03692 


0.81896 


0.32070 


0.07079 


0.81896 


0.32289 


0.04124 


0.81896 


0.33535 


0.07111 


0.81896 


0.33851 


0.04599 


0.81896 


0.34972 


0.07177 


0.81896 


0.34740 


0.04889 


0.818% 


0.35778 


0.07224 


0.81896 


0.35399 


0.05110 


0.81896 


0.36374 


0.07270 


0.81896 


0.36933 


0.05667 


0.81896 


0.37744 


0.07400 


0.81896 


0.38450 


0.06262 


0.81896 


0.39081 


0.07568 


0.81896 


0.38643 


0.06344 


0.81896 


0.39248 


0.07591 


0.81896 


0.39950 


0.06905 


0.81896 


0.40385 


0.07778 


0.81896 


0.41431 


0.07593 


0.81896 


0.41655 


0.08030 


0.81896 


0.42889 


0.08323 


0.81896 


0.42889 


0.08323 


35 






35 






0.84998 


0.00000 


0.00001 


0.84998 


0.00000 


0.10925 


0.84998 


0.01629 


0.00084 


0.84998 


0.01870 


0.10925 


0.84998 


0.03258 


0.00171 


0.84998 


0.03735 


0.10925 
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0.84998 


0.04887 


0.00262 


0.84998 


0.05595 


0.10925 


0.84998 


0.06515 


0.00362 


0.84998 


0.07450 


0.10925 


0.84998 


0.06824 


0.00382 


0.84998 


0.07801 


0.10925 


0.84998 


0.06839 


0.00383 


0.84998 


0.07818 


0.10925 


0.84998 


0.06854 


0.00384 


0.84998 


0.07828 


0.08202 


0.84998 


0.08143 


0.00470 


0.84998 


0.08995 


0.08180 


0.84998 


0.09770 


0.00582 


0.84998 


0.10769 


0.08133 


0.84998 


0.11397 


0.00711 


0.84998 


0.12530 


0.08070 


0.84998 


0.13022 


0.00848 


0.84998 


0.14279 


0.08001 


0.84998 


0.14647 


0.01000 


0.84998 


0.16013 


0.07928 


0.84998 


0.16269 


0.01165 


0.84998 


0.17729 


0.07846 


0.84998 


0.17890 


0.01350 


0.84998 


0.19429 


0.07767 


0.84998 


0.19509 


0.01553 


0.84998 


0.21106 


0.07680 


0.84998 


0.21125 


0.01777 


0.84998 


0.22764 


0.07600 


0.84998 


0.22737 


0.02024 


0.84998 


0.24398 


0.07527 


0.84998 


0.24345 


0.02297 


0.84998 


0.26008 


0.07460 


0.84998 


0.25948 


0.02598 


0.84998 


0.27592 


0.07406 


0.84998 


0.27377 


0.02893 


0.84998 


0.28986 


0.07371 


0.84998 


0.27546 


0.02930 


0.84998 


0.29150 


0.07368 


0.84998 


0.29135 


0.03291 


0.84998 


0.30679 


0.07348 


0.84998 


0.30717 


0.03690 


0.84998 


0.32179 


0.07352 


0.84998 


0.32290 


0.04125 


0.84998 


0.33648 


0.07381 


0.84998 


0.33851 


0.04599 


0.84998 


0.35085 


0.07439 


0.84998 


0.34740 


0.04888 


0.84998 


0.35895 


0.07486 


0.84998 


0.35399 


0.05110 


0.84998 


0.36491 


. 0.07529 


0.84998 


0.36934 


0.05667 


0.84998 


0.37863 


0.07654 


0.84998 


0.38451 


0.06264 


0.84998 


0.39202 


0.07817 


0.84998 


0.38643 


0.06344 


0.84998 


0.39369 


0.07840 


0.84998 


0.39950 


0.06904 


0.84998 


0.40506 


0.08020 


0.84998 


0.41431 


0.07593 


0.84998 


0.41776 


0.08266 


0.84998 


0.42889 


0.08323 


0.84998 


0.43012 


0.08556 


0.84998 


0.43654 


0.08728 


0.84998 


0.43654 


0.08728 


36 






36 






0.87780 


0.00000 


0.00001 


0.87780 


0.00000 


0.10925 


0.87780 


0.01629 


0.00083 


0.87780 


0.01870 


0.10925 


0.87780 


0.03258 


0.00169 


0.87780 


0.03735 


0.10925 


0.87780 


0.04887 


0.00263 


0.87780 


0.05595 


0.10925 


0.87780 


0.06515 


0.00363 


0.87780 


0.07450 


0.10925 


0.87780 


0.06824 


0.00382 


0.87780 


0.07801 


0.10925 


0.87780 


0.06839 


0.00383 


0.87780 


0.07818 


0.10925 


0.87780 


0.06854 


0.00384 


0.87780 


0.07828 


0.07371 


0.87780 


0.08143 


0.00468 


0.87780 


0.08906 


0.07377 


0.87780 


0.09770 


0.00585 


0.87780 


0.10668 


0.07368 


0.87780 


0.11397 


0.00711 


0.87780 


0.12420 


0.07350 


0.87780 


0.13022 


0.00845 


0.87780 


0.14161 


0.07327 


0.87780 


0.14647 


0.00999 


0.87780 


0.15887 


0.07292 


0.87780 


0.16270 


0.01166 


0.87780 


0.17600 


0.07253 


0.87780 


0.17890 


0.01348 


0.87780 


0.19298 


0.07221 
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0.87780 


0.19509 


0.01553 


0.87780 


0.20974 


0.07173 


0.87780 


0.21125 


0.01778 


0.87780 


0.22633 


0.07134 


0.87780 


0.22737 


0.02025 


0.87780 


0.24271 


0.07103 


0.87780 


0.24346 


0.02298 


0.87780 


0.25886 


0.07081 


0.87780 


0.25949 


0.02599 


0.87780 


0.27478 


0.07070 


0.87780 


0.27377 


0.02894 


0.87780 


0.28880 


0.07076 


0.87780 


0.27546 


0.02931 


0.87780 


0.29044 


0.07075 


0.87780 


0.29136 


0.03294 


0.87780 


0.30584 


0.07098 


0.87780 


0.30717 


0.03691 


0.87780 


0.32095 


0.07142 


0.87780 


0.32290 


0.04126 


0.87780 


0.33573 


0.07202 


0.87780 


0.33851 


0.04598 


0.87780 


0.35028 


0.07308 


0.87780 


0.34740 


0.04889 


0.87780 


0.35844 


0.07373 


0.87780 


0.35399 


0.05112 


0.87780 


0.36448 


0.07433 


0.87780 


0.36934 


0.05668 


0.87780 


0.37837 


0.07599 


0.87780 


0.38451 


0.06265 


0.87780 


0.39188 


0.07789 


0.87780 


0.38643 


0.06344 


0.87780 


0.39356 


0.07814 


0.87780 


0.39950 


0.06904 


0.87780 


0.40606 


0.08221 


0.87780 


0.41431 


0.07593 


0.87780 


0.41863 


0.08436 


0.87780 


0.42890 


0.08323 


0.87780 


0.43089 


0.08703 


0.87780 


0.43654 


0.08729 


0.87780 


0.43728 


0.08867 


0.87780 


0.44325 


0.09098 


0.87780 


0.44325 


0.09098 
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36 






0.93164 


0.00000 


0.00002 


0.93164 


0.00000 


0.10925 


0.93164 


0.01629 


0.00084 


0.93164 


0.01870 


0.10925 


0.93164 


0.03258 


0.00171 


0.93164 


0.03735 


0.10925 


0.93164 


0.04887 


0.00263 


0.93164 


0.05595 


0.10925 


0.93164 


0.06515 


0.00363 


0.93164 


0.07450 


0.10925 


0.93164 


0.06824 


0.00383 


0.93164 


0.07801 


0.10925 


0.93164 


0.06839 


0.00384 


0.93164 


0.07818 


0.10925 


0.93164 


0.06854 


0.00385 


0.93164 


0.07828 


0.05684 


0.93164 


0.08143 


0.00469 


0.93164 


0.08723 


0.05716 


0.93164 


0.09770 


0.00583 


0.93164 


0.10457 


0.05778 


0.93164 


0.11397 


0.00711 


0.93164 


0.12184 


0.05822 


0.93164 


0.13022 


0.00847 


0.93164 


0.13906 


0.05875 


0.93164 


0.14647 


0.01000 


0.93164 


0.15616 


0.05915 


0.93164 


0.16269 


0.01165 


0.93164 


0.17319 


0.05970 


0.93164 


0.17890 


0.01350 


0.93164 


0.19008 


0.06011 


0.93164 


0.19509 


0.01553 


0.93164 


0.20685 


0.06063 


0.93164 


0.21125 


0.01777 


0.93164 


0.22350 


0.06129 


0.93164 


0.22737 


0.02024 


0.93164 


0.23998 


0.06199 


0.93164 


0.24345 


0.02298 


0.93164 


0.25628 


0.06281 


0.93164 


0.25949 


0.02599 


0.93164 


0.27242 


0.06382 


0.93164 


0.27377 


0.02893 


0.93164 


0.28666 


0.06482 


0.93164 


0.27546 


0.02930 


0.93164 


0.28835 


0.06496 


0.93164 


0.29136 


0.03293 


0.93164 


0.30406 


0.06632 


0.93164 


0.30718 


0.03691 


0.93164 


0.31957 


0.06796 


0.93164 


0.32291 


0.04127 


0.93164 


0.33484 


0.06987 


0.93164 


0.33852 


0.04600 


0.93164 


0.34986 


0.07209 
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0.93164 


0.34740 


0.04887 


0.93164 


0.35835 


0.07352 


0.93164 


0.35400 


0.05113 


0.93164 


0.36458 


0.07456 


0.93164 


0.36933 


0.05666 


0.93164 


0.37909 


0.07753 


0.93164 


0.38451 


0.06265 


0.93164 


0.39323 


0.08068 


0.93164 


0.38642 


0.06341 


0.93164 


0.39505 


0.08121 


0.93164 


0.39951 


0.06906 


0.93164 


0.40798 


0.08607 


0.93164 


0.41431 


0.07593 


0.93164 


0.42031 


0.08764 


0.93164 


0.42890 


0.08323 


0.93164 


0.43239 


0.08988 


0.93164 


0.43654 


0.08729 


0.93164 


0.43870 


0.09135 


0.93164 


0.44325 


0.09098 


0.93164 


0.44446 


0.09322 


0.93816 


0.45582 


0.09828 


0.93164 


0.45582 


0.09828 


38 






38 






1.00000 


0.00000 


0.00001 


1.00000 


0.00000 


0.03000 


1.00000 


0.01629 


0.00083 


1.00000 


0.01694 


0.03000 


1.00000 


0.03258 


0.00170 


1.00000 


0.03384 


0.03000 


1.00000 


0.04887 


0.00263 


1.00000 


0.05069 


0.03000 


1.00000 


0.06515 


0.00362 


1.00000 


0.06749 


0.03000 


1.00000 


0.06824 


0.00381 


1.00000 


0.07067 


0.03000 


1.00000 


0.06839 


0.00382 


1.00000 


0.07116 


0.03371 


1.00000 


0.06854 


0.00383 


1.00000 


0.07132 


0.03372 


1.00000 


0.08143 


0.00468 


1.00000 


0.08473 


0.03452 


1.00000 


0.09770 


0.00584 


1.00000 


0.10164 


0.03560 


1.00000 


0.11397 


0.00710 


1.00000 


0.11854 


0.03676 


1.00000 


0.13022 


0.00847 


1.00000 


0.13542 


0.03803 


1.00000 


0.14647 


0.00999 


1.00000 


0.15227 


0.03943 


1.00000 


0.16269 


0.01165 


1.00000 


0.16910 


0.04097 


1.00000 


0.17890 


0.01349 


1.00000 


0.18590 


0.04268 


1.00000 


0.19509 


0.01552 


1.00000 


0.20266 


0.04457 


1.00000 


0.21125 


0.01777 


1.00000 


0.21938 


0.04665 


1.00000 


0.22737 


0.02024 


1.00000 


0.23605 


0.04897 


1.00000 


0.24345 


0.02298 


1.00000 


0.25265 


0.05154 


1.00000 


0.25949 


0.02599 


1.00000 


0.26920 


0.05438 


1.00000 


0.27377 


0.02894 


1.00000 


0.28392 


0.05717 


1.00000 


0.27546 


0.02930 


1.00000 


' 0.28566 


0.05752 


1.00000 


0.29136 


0.03294 


1.00000 


0.30203 


0.06098 


1.00000 


0.30718 


0.03692 


1.00000 


0.31830 


0.06479 


1.00000 


0.32290 


0.04127 


1.00000 


0.33445 


0.06895 


1.00000 


0.33851 


0.04600 


1.00000 


0.35047 


0.07351 


1.00000 


0.34740 


0.04889 


1.00000 


0.35959 


0.07630 


1.00000 


0.35400 


0.05113 


* 1.00000 


0.36634 


0.07847 


1.00000 


0.36934 


0.05668 


1.00000 


0.38205 


0.08384 


1.00000 


0.38452 


0.06265 


1.00000 


0.39757 


. 0.08966 


1.00000 


0.38643 


0.06344 


1.00000 


0.39952 


0.09042 


1.00000 


0.39951 


0.06907 


• 1.00000 


0.41045 


0.09102 


1.00000 


0.41431 


0.07593 


1.00000 


0.42245 


0.09180 


1.00000 


0.42890 


0.08323 


1.00000 


0.43429 


0.09349 


1.00000 


0.43654 


0.08729 


1.00000 


0.44052 


0.09475 


1.00000 


0.44325 


0.09099 


1.00000 


0.44599 


0.09608 
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1.00000 


0.45736 


0.09918 


1.00000 


0.45756 


0.09954 


1.00000 


0.47120 


0.10780 


1.00000 


0.47120 


0.10780 


14 






14 






1.01500 


0.00000 


0.01500 


1.01500 


0.00000 


0.01500 


1.01500 


0.01694 


0.01541 


1.01500 


0.01694 


0.01541 


1.01500 


0.03258 


0.01585 


1.01500 


0.03258 


0.01585 


1.01500 


0.07000 


0.01692 


1.01500 


0.07000 


0.01692 


1.01500 


0.10000 


0.02070 


1.01500 


0.10000 


0.02070 


1.01500 


0.13250 


0.02330 


1.01500 


0.13250 


0.02330 


1.01500 


0.18500 


0.02859 


1.01500 


0.18500 


0.02859 


1.01500 


0.25500 


0.03876 


1.01500 


0.25500 


0.03876 


1.01500 


0.30400 


0.04895 


1.01500 


0.30400 


0.04895 


1.01500 


0.32290 


0.05407 


1.01500 


0.32290 


0.05407 


1.01500 


0.38300 


0.07324 


1.01500 


0.38300 


0.07324 


1.01500 


0.42600 


0.08751 


1.01500 


0.42600 


0.08751 


1.01500 


0.44450 


0.09353 


1.01500 


0.44450 


0.09353 


1.01500 


0.47120 


0.10780 


1.01500 


0.47120 


0.10780 
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APPENDIX B. SAMPLE INPUT FILES 



Mach 3.0 (Inviscid) Input file. Start 



$ GLOBAL 

CHIMRA= .F., NSTEPS=200, RESTRT= . F. , NSAVE =10, 
NQT = 0, 

$END 

$FLOINP 

ALPHA = 0, FSMACH = 0.3, REY = 3.19E6, TINF = 
520.000, 

$END 

$VARGAM 

$END 

' $GRDNAM 

NAME = 'Waverider Half-Body, 42x45x50 grid'. 



$END 

$NITERS 

$END 

$METPRM 



IRHS = 0, 
$END 
$TIMACU 


ILHS = 2, 


IDISS 


= 2, 


DT = 0. 

$END 
$SMOACU 


01, ITIME= 1, 


TFOSO = 


= 1.00, 


ISPEC J= 2, 


DIS2J = 0.0, 


DIS4J 


C\1 

O 

II 


ISPECK= 2, 


DIS2K = 0.0, 


DIS4K 


K. 

CM 

O 

II 


ISPECL= 2, 


DIS2L = 0.0, 


DIS4L 


CM 

O 

II 



CFLMIN=0. 0, 



SMOO 

EPSE 

$END 

$VISINP 

VISCJ 

NTURB 

ITTYP 

ITDIR 

JTLS 

JTLE 

KTLS 

KTLE 

LTLS 

LTLE 



1 . 00 , 

0.35, 



F., 

0 , 

1 , 

3, 

1 , 

42, 

1, 

45, 

1 , 

50, 



VISCK = . F. , VISCL = . F. 



TLPAR1= 0.3, 
$END 
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$BCINP 



NBC 


= 


6 , 




IBTYP 


= 


1, 


14, 


IBDIR 


= 


3, 


2, 


JBCS 


= 


1, 


1, 


JBCE 


= 


42, 


42, 


KBCS 


= 


1, 


1, 


KBCE 


= 


45, 


1, 


LBCS 


= 


1, 


1, 


LBCE 

$END 


= 


1, 


49, 



$SCEINP 

$END 



12, 


54, 


12, 


32 


-1, 


-2, 


1, 


-3 


42, 


1, 


1, 


1 


42, 


42, 


1, 


42 


1, 


45, 


1, 


1 


45, 


45, 


45, 


45 


1, 


1, 


1, 


50 


50, 


50, 


50, 


50 
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Mach 3.0 (Viscous) Input File, Restart 



$ GLOBAL 

CHIMRA= 
=100, NQT 
$END 
$FLOINP 

ALPHA = 
520.000, 

$END 

$VARGAM 

$END 

$GRDNAM 

NAME = 

$END 

$NITERS 

$END 

$METPRM 

IRHS = 
$END 
$TIMACU 
DT 

$END 

$SMOACU 

ISPECJ= 
ISPECK= 
ISPECL= 
SMOO = 
EPSE = 
$END 
$VISINP 

VISCJ = 
NTURB = 
ITTYP = 
ITDIR = 
JTLS = 
JTLE = 
KTLS = 
KTLE = 
LTLS = 
LTLE = 
TLPAR1= 
$END 
$BCINP 
NBC 



. F., NSTEPS=10000, RESTRT= .T., NSAVE 

0, 

0, FSMACH = 0.3, REY = 3.19E6, TINF = 



Waverider Half-Body, 4.2x45x50 grid'. 



0, ILHS = 2, IDISS = 2, 



0.07, ITIME= 1, TFOSO = 1.00, CFLMIN=0.01, 



2, DIS2J = 0.0, DIS4J 
2, DIS2K = 0.0, DIS4K 
2, DIS2L = 0.0, DIS4L 
1 . 00 , 

0.35, 



0 . 2 , 

0 . 2 , 

0 . 2 , 



.T., VTSCK = .F., VISCL = -T., 

1 , 

1 , 

3, 

1 , 

42, 

1 , 

45, 

1 , 

50, 

0.3, 



6 , 
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IBTYP 


= 


5, 


14, 


12, 


54, 


12, 


32 


IBDIR 


= 


3, 


2, 


-1, 


-2, 


1, 


-3 


JBCS 


= 


1, 


1, 


42, 


1, 


1, 


1 


JBCE 


= 


42, 


42, 


42, 


42, 


1, 


42 


KBCS 


= 


1, 


1, 


1, 


45, 


1, 


1 


KBCE 


= 


45, 


1, 


45, 


45, 


45, 


45 


LBCS 


— 


1, 


1, 


1, 


1, 


1, 


50 


LBCE 


= 


1/ 


49, 


50, 


50, 


50, 


50 



$END 

$SCEINP 

$END 
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Mach 6.0 (Inviscid) Input file. Start 

$ GLOBAL 

CHIMRA= . F., NSTEPS=100 , RESTRT= . F. , NSAVE =10, 
NQT = 0, 

$END 
$ FLOINP 

ALPHA = 0, FSMACH = 6.0, REY = 6.39E7, TINF = 
520.000, 

$END 

$VARGAM 

$END 

$GRDNAM 

NAME = 'Waverider Half-Body, 42x45x50 grid’. 



$END 

$NITERS 

$END 

$METPRM 



IRHS = 0, 
$END 


ILHS = 2, 


IDISS 


= 2, 


$TIMACU 








DT = 0. 

$END 


001, ITIME= 1, 


TFOSO 


= 1.00 


$SMOACU 








ISPECJ= 2, 


DIS2J = 2.0, 


DIS4J 


= 0.2, 


ISPECK= 2, 


DIS2K = 2.0, 


DIS4K 


CM 

O 

II 


ISPECL= 2, 


DIS2L = 2.0, 


DIS4L 


= 0.2, 



CFLMIN=0. 0, 



SMOO 
EPSE 
$END 
5VISINP 
VISCJ 
NTURB 
ITTYP 
ITDIR 
JTLS 
JTLE 
KTLS 
KTLE 
LTLS 
LTLE 
TLPAR1= 0.3, 
$END 
$BCINP 

NBC = 6, 

IBTYP = 



1 . 00 , 

0.35, 



F., 

0 , 

1 , 

3, 

1 , 

42, 

1 , 

45, 

1 , 

50, 



VISCK = .F., VISCL = .F., 



1, 14, 



12, 54, 12, 32, 
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IBDIR = 


3, 


2, 


JBCS = 


1, 


1, 


JBCE = 


42, 


42, 


KBCS = 


1, 


1, 


KBCE = 


45, 


1, 


LBCS = 


1, 


1, 


LBCE = 
$END 
$SCEINP 


1, 


50, 



$END 



-1, 


-2, 


1, 


-3 


42, 


1, 


1, 


1 


42, 


42, 


1, 


42 


1, 


45, 


1, 


1 


45, 


45, 


45, 


45 


1, 


1, 


1, 


50 


50, 


50, 


50, 


50 
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Mach 6.0 (Viscous) File input, Restart 



• F., 
0, 



NSTEPS=4000, RESTRT= .T., 



NSAVE 



$GLOBAL 

CHIMRA^ 

= 100 , 

NQT 
$END 
$FLOINP 

ALPHA = 0, FSMACH = 6.0, REY = 6.39E7, TINF = 
520.000, 

$END 

$VARGAM 

$END 

$GRDNAM 

NAME = 'Waverider Half-Body, 42x45x50 grid', 

$END 

$NITERS 

$END 

$METPRM 



01 , 



IRHS = 0, 


ILHS = 2, 


IDISS = 


= 2, 


$END 

$TIMACU 

DT = 0. 


005, ITIME= 1, 


TFOSO = 


1.00, CFLMIN : 


$END 

$SMOACU 

ISPECJ= 2, 


DIS2J = 2.0, 


DIS4J = 


0.2, 


ISPECK= 2, 


DIS2K = 2.0, 


DIS4K = 


0.2, 


ISPECL= 2, 


DIS2L =2.0, 


DIS4L = 


0.2, 


SMOO = 1 . 
EPSE = 0. 
$END 
$VISINP 

VISCJ = .T 


00, 

35, 

. , VISCK = . F. 


, VISCL 


= .T., 



NTURB 

ITTYP 

ITDIR 

JTLS 

JTLE 

KTLS 

KTLE 

LTLS 

LTLE 



1 , 

1 , 

3, 

1 , 

42, 

1 , 

45, 

1 , 

50, 



TLPAR1= 0.3, 
$END 
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$BCINP 



NBC 


= 


6, 




IBTYP 


= 


5, 


14, 


IBDIR 


= 


3, 


2, 


JBCS 


= 


1, 


If 


JBCE 


= 


42, 


42, 


KBCS 


= 


lr 


If 


KBCE 


= 


45, 


If 


LBCS 


= 


lr 


If 


LBCE 

$END 

$SCEINP 




1, 


O 

LO 



$END 



12, 


54, 


12, 


32 


-If 


-2, 


If 


-3 


42, 


1, 


If 


1 


42, 


42, 


If 


42 


If 


45, 


If 


1 


45, 


45, 


45, 


45 


If 


1, 


1, 


50 


50, 


50, 


50, 


50 
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Mach 6.0 (Inviscid) Engine Flow-Through model Input file, Start 



$ GLOBAL 

CHIMRA= . F., NSTEPS=100, RESTRT= .F., NSAVE =10, 
NQT = 0, 

$END 

$FLOINP 

ALPHA = 0, FSMACH = 6.0, REY = 6.39E7, TINF = 
520.000, 

$END 

$VARGAM 

$END 

$GRDNAM 

NAME = 'Waverider Half-Body, 42x45x50 grid', 



$END 

$NITERS 

$END 

$METPRM 

IRHS = 
$END 
$TIMACU 
DT 

$END 

$SMOACU 

ISPEGJ= 
ISPECK= 
ISPECL= 
SMOO = 
EPSE = 
$END 
$VISINP 

VISCJ = 
NTURB = 
ITTYP = 
ITDIR = 
JTLS = 
JTLE = 
KTLS = 
KTLE = 
LTLS = 
LTLE = 
TLPAR1= 
$END 
$BCINP 



0, ILHS = 2, 



0.001, ITIME= 1, 



2, DIS2J = 2.0, 
2, DIS2K = 2.0, 
2, DIS2L = 2.0, 
1 . 00 , 

0.35, 



. F. , VISCK = . F 

0 , 

1 , 

3, 

1 , 

42, 

1 , 

45, 

1 , 

50, 

0.3, 



IDISS = 2, 

TFOSO = 1.00, 

DIS4J = 0.2, 
DIS4K = 0.2, 
DIS4L = 0.2, 

, VISCL = . F. 



CFLMIN=0. 0, 
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8 



NBC 



IBTYP 


= 


1, 


\ — i 


12, 


54, 


12, 


32, 


30, 


40, 


IBDIR 


= 


3, 


2, 


-1, 


-2, 


1, 


-3, 


3, 


3, 


JBCS 


= 


1, 


1, 


42, 


1, 


1, 


1, 


38, 


38, 


JBCE 


= 


42, 


42, 


42, 


42, 


1, 


42, 


42, 


42, 


KBCS 


= 


1, 


1, 


1, 


45, 


1, 


1, 


21, 


37, 


KBCE 


= 


45, 


1, 


45, 


45, 


45, 


45, 


26, 


40, 


LBCS 


= 


1, 


1, 


1, 


1, 


1, 


50, 


1, 


1, 


LBCE 


= 


1, 


50, 


50, 


50, 


50, 


o 

LO 


1, 


1, 



$END 

$SCEINP 

$END 
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APPENDIX C. FLOW SOLUTIONS 




Figure C-l. Mach 0.3 Inviscid Convergence History 
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Residual Residual 




Figure C-3. Mach 0.4 Inviscid Convergence History 




Figure C-4. Mach 0.4 Viscous Convergence History' 
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Residual Res dual 




Figure C-5. Mach 0.5 Inviscid Convergence History 




Figure C-6. Mach 0.5 Viscous Convergence History 
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Figure C-7. Mach 0.6 Inviscid Convergence History 




Figure C-8. Mach 0.6 Viscous Convergence History 
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Residual Ftesdual 




Figure C-9. Mach 0.7 Inviscid Convergence History 




Figure C-10. Mach 0.7 Viscous Convergence History 
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Residual Residual 




Figure C-l 1. Mach 0.8 Inviscid Convergence History 




Figure C-12. Mach 0.8 Viscous Convergence History 
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Figure C-13. Mach 0.9 Inviscid Convergence History 




Figure C-14. Mach 0.9 Viscous Convergence History 
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Residual Residual 




Figure C-15. Mach 1.2 Inviscid Convergence History 




Figure C-16. Mach 1.2 Viscous Convergence History 
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Figure C-17. Mach 2.0 Inviscid Convergence History 




Figure C-18. Mach 2.0 Viscous Convergence History 
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Residual Residual 




Figure C-19. Mach 4.0 Inviscid Convergence History 




Figure C-20. Mach 4.0 Viscous Convergence History 
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FtosidLBl Residual 




Figure C-21. Mach 6.0 Inviscid Convergence History 




Figure C-22. Mach 6.0 Viscous Convergence History 
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Figure C-23. Mach 6.0 Inviscid Engine Flow-Through Convergence History 
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Figure C-24. Force Coefficients vs. Iterations Mach 0.3 Inviscid 
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Figure C-25. Force Coefficients vs. Iterations Mach 0.3 Viscous 
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Figure C-26. Force Coefficients vs. Iterations Mach 0.4 Inviscid 
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Figure C-27. Force Coefficients vs. Iterations Mach 0.4 Viscous 
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Figure C-28. Force Coefficients vs. Iterations Mach 0.5 Inviscid 
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Figure C-29. Force Coefficients vs. Iterations Mach 0.5 Viscous 
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Figure C-30. Force Coefficients vs. Iterations Mach 0.6 Inviscid 
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Figure C-31. Force Coefficients vs. Iterations Mach 0.6 Viscous 



102 





Figure C-32. Force Coefficients vs. Iterations Mach 0.7 Inviscid 
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Figure C-33. Force Coefficients vs. Iterations Mach 0.7 Viscous 
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Figure C-34. Force Coefficients vs. Iterations Mach 0.8 Inviscid 



105 





Figure C-35. Force Coefficients vs. Iterations Mach 0.8 Viscous 
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Figure C-36. Force Coefficients vs. Iterations Mach 0.9 Inviscid 
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Figure C-37. Force CoefTicients vs. Iterations Mach 0.9 Viscous 
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Figure C-38. Force Coefficients vs. Iterations Mach 1.2 Inviscid 
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Figure C-41. Force Coefficients vs. Iterations mach 2.0 Viscous 
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Figure C-42. Force Coefficients vs. Iterations Mach 4.0 Inviscid 
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Figure C-43. Force Coefficients vs. Iterations Mach 4.0 Viscous 
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Figure C-44. Force Coefficients vs. Iterations Mach 6.0 Inviscid 
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Figure C-45. Force Coefficients vs. Iterations Mach 6.0 Viscous 
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Figure C-46. Force Coefficients vs. Iterations Mach 6.0 Engine Flow-Through 
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Figure C-47. Mach 0.3 Inviscid Normalized Surface Pressure and Mach Contours 




Figure C-48. Mach 0.3 Viscous Normalized Surface Pressure and Mach Contours 
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Figure C-50. Mach 0.4 Viscous Normalized Surface Pressure and Mach Contours 
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Figure C-51. Mach 0.5 Inviscid Normalized Surface Pressure and Mach Contours 




Figure C-52. Mach 0.5 Viscous Normalized Surface Pressure and Mach Contours 
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Figure C-53. Mach 0.6 Inviscid Normalized Surface Pressure and Mach Contours 




Figure C-54. Mach 0.6 Viscous Normalized Surface Pressure and Mach Contours 
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Figure C-55. Mach 0.7 Inviscid Normalized Surface Pressure and Mach Contours 




Figure C-56. Mach 0.7 Viscous Normalized Surface Pressure and Mach Contours 
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Figure C-58. Mach 0.8 Viscous Normalized Surface Pressure and Mach Contours 
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Figure C-59. Mach 0.9 Inviscid Normalized Surface Pressure and Mach Contours 




0.423 1.028 1.633 0.00 0.73 1.45 



Figure C-60. Mach 0.9 Viscous Normalized Surface Pressure and Mach Contours 
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Figure C-61. Mach 1.2 Inviscid Normalized Pressure and Mach Contours 




Figure C-62. Mach 1.2 Viscous Normalized Surface Pressure and Mach Contours 
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Figure C-63. Mach 2.0 Inviscid Normalized Surface Pressure and Mach Contours 
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Figure C-66. Mach 4.0 Viscous Normalized Surface Pressure and Mach Contours 
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Figure C-67. Mach 6.0 Inviscid Normalized Surface Pressure and Mach Contours 




Figure C-68. Mach 6.0 Viscous Normalized Surface Pressure and Mach Contours 
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Figure C-69. Mach 6.0 Engine Flow-Through Normalized Surface Pressure and Mach 

Contours 
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